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X1908+075: A PULSAR ORBITING IN THE STELLAR WIND OF A MASSIVE COMPANION 
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ABSTRACT 

We have observed the persistent but optically unidentified X-ray source X1908+075 with the 
PC A and HEXTE instruments on RXTE. The binary nature of this source was established by 
IWen. Remillard. fc Bradl lj2000fl who found a 4.4-day orbital period in results from the RXTE ASM. 
We report the discovery of 605 s pulsations in the X-ray flux. The Doppler delay curve is measured 
and provides a mass function of 6.1 Mq which is a lower limit to the mass of the binary companion 
of the neutron star. The degree of attenuation of the low-energy end of the spectrum is found to be a 
strong function of orbital phase. A simple model of absorption in a stellar wind from the companion 
star fits the orbital phase dependence reasonably well and limits the orbital inclination angle to the 
range 38° — 72°. These measured parameters lead to an orbital separation of ~ 60 — 80 lt-s, a mass for 
the companion star in the range 9-31 Mq, and an upper limit to the size of the companion of ~ 22 Rq. 
From our analysis we also infer a wind mass loss rate from the companion star of > 1.3 x 10 -6 Mq 
yr -1 and, when the properties of the companion star and the effects of photoionization are considered, 
likely > 4 x 10~ 6 Mq yr _1 . Such a high rate is inconsistent with the allowed masses and radii that we 
find for a main sequence or modestly evolved star unless the mass loss rate is enhanced in the binary 
system relative to that of an isolated star. We discuss the possibility that the companion might be a 
Wolf-Rayet star that could evolve to become a black hole in 10 4 to 10 5 yr. If so, this would be the 
first identified progenitor of a neutron star-black hole binary. 

Subject headings: X-rays: binaries — pulsars: general — pulsars: individual (X1908+075) — stars: 
evolution — stars: winds, outflows 



1. INTRODUCTION 

X1908+075 is an optically unidentified, highly ab- 
sorbed, and relatively faint X-ray source that appeared 
in surveys carried out with instruments on the Uhuru, 
OSO 7, Ariel 5, HEAO-1, and EXOSAT satellites. The 
early detect ions and position determi n ations are sum- 
marized by IWen. Remillard. fc Bradtl (|2000, hereafter 
WRB). They conclude that the position of X1908+075 
is likely to be within the overlapping region of the error 
box of a source detected in an Einstein IPC image and 
one of an array of HE A O 1 A-3 position "diamonds" , and 
is thus known with an accuracy of ~ 1'. Inspection of 
the POSS plates within the source error box reveals no 
optical counterpart down to magnitude 20. This is con- 
sistent with the heavy optical extinction implied by the 
interstellar hydrogen column density of ~ 4 x 10 22 atoms 
cm~ 2 measured on the basi s of the low-energy absorptio n 
in the X-ray spectrum fsee lShull fc van Steenber J l985L 

The intensity of X1908+075 has been monitored for 
the past eight years using the All-Sky Monitor (ASM) 
aboard the Rossi X-Ray Timing Explorer (RXTE), and 
has typically been in the range 2-12 mCrab in the 2-12 
keV energy band. WRB analyzed data from the first 
three years of operation of the ASM, and thereby dis- 
covered a 4.4 day periodicity in the X-ray intensity. The 
periodic component of the intensity variations is clearly 
energy dependent in both strength and detailed form. At 
most orbital phases, the variation is roughly sinusoidal 
while a relatively sharp dip forms the minimum in the 
5-12 keV band. These characteristics suggest that the 
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modulation is produced by a varying amount of absorp- 
tion along the line of sight as the source moves through 
the stellar wind of a massive companion star. The hard 
X-ray spectrum led WRB to also suggest that this source 
could be an X-ray pulsar. 

The possibility that X1908+075 might be an X-ray 
pulsar led us to carry out a set of pointed observations 
with the PCA and HEXTE instruments on RXTE in late 
2000 and early 2001. The data revealed the presence of 
strong X-ray pulsations at a period of 605 seconds. We 
were also able to detect Doppler delays in the pulse ar- 
rival times. However, because the number of indepen- 
dent high-quality pulse arrival times obtained from this 
data set was small, it proved difficult to unambiguously 
disentangle orbital effects from intrinsic changes in the 
pulse period. The latter could, in principle, be quite 
large for a neutron star rotating with a period as long as 
600 s (see, e.g.. lBildsten et al.ll997tlDelgado-Marti et alJ 
2001). Therefore, we obtained additional observations of 
the source with RXTE during late 2002 and early 2003. 

In this paper we report the results of our analysis of 
the RXTE observations of X1908+075. The pointed ob- 
servations are described in §2. A pulse timing analysis is 
described in §3. We present an accurate pulse period and 
a determination of the spin-down of the neutron star due 
to accretion and magnetic torques. The orbital Doppler 
delay curve is measured, thereby confirming the 4.4-day 
period found in the ASM X-ray light curve. The resul- 
tant mass function is 6.1 Mq, indicating that the pulsar 
does indeed orbit a massive companion star. The results 
of an orbital-phase-dependent spectral analysis are pre- 
sented in §4. We detect a very pronounced modulation 
of the low energy attenuation as a function of orbital 
phase. In §5 we model this modulation by absorption in 
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a spherically symmetric stellar wind, whereby we obtain 
constraints on the orbital inclination, on the properties 
of the companion star, and on the stellar wind. We dis- 
cuss the implications of our results in §6, including the 
possibility that this system may be the progenitor of a 
neutron star-black hole binary. 

2. OBSERVATIONS 

Pointed observations of X1908+075 were made with 
the Proportional Counter Array (PCA) and High-Energy 
X-ray Timing Experiment (HEXTE) on RXTE. The 
PCA consists of 5 Proportional Counter Units (PCUs) 
that are sensitive to X-ray photons in the range 2.5- 
60 keV. Each PCU has a collecting area of ~ 1400 cm 2 
and a collimator to limit the field of view to 1° in radius 
l)Jahoda et al.lll996T) . Some of the PCUs are operated 
with reduced duty cycle in order to avoid problems asso- 
ciated with constant use. Only PCUs and 2 were used 
in every observation. Data from the PCA observations 
were telemetered to the ground in the "GoodXenon" 
mode, which includes information on each good event 
with 1 /is time resolution and the instrument's full en- 
ergy resolution (255 channels). The HEXTE comprises 
two clusters, each of which includes 4 Nal scintillation 
detectors sensitive to photons in the range 15-250 keV 
that provide a total collecting area of 800 cm 2 per clus- 
ter. The d etectors in each cluster view a common 1° 
radius field l|Rothschild et alJll998|) . 

Between 2000 November 19 and 2001 February 22, we 
obtained 15 short observations with a typical exposure 
time in each of ~ 2000 s. A total exposure of 35 ks was 
obtained in these observations, which we refer to col- 
lectively as "epoch 1" . Over all of the 15 observations 
there were, on average, 3.0 PCUs in operation. Eight 
of the observations were made during one orbital cycle 
of X1908+075 (~ 4 days), while the remaining 7 were 
performed nearly 3 months later. Fourier transforms of 
these data clearly exhibited the presence of X-ray pulsa- 
tions with a period of 605 s. 

Additional observations with the PCA and HEXTE 
were carried out on 39 occasions from 2002 December 23 
through 2003 January 5 and on 32 occasions from 2003 
January 30 through 2003 February 8 ("epoch 2"; 196 ks 
total exposure). On average, 3.5 PCUs were in opera- 
tion. Each of these 71 observations in epoch 2 typically 
yielded ~ 3000 s of net exposure most often in the form 
of two contiguous time intervals separated by a ~ 2000 s 
gap due to Earth occultation of the target. Raw count- 
ing rate data from two of these 71 pointed observations 
are shown in Fig. ^to illustrate the typical appearance 
of the X-ray pulsations and accompanying source vari- 
ability. The X-ray intensity, in addition to exhibiting 
obvious pulsations, is also quite variable on timescales 
comparable to the pulse period and longer. 

For pulsar timing purposes we used events with ener- 
gies in the range 3.7-17 keV (channels 3-40) from both 
left and right sections of layer 1 of all operating PCUs. 
Event times were reduced to the Solar System barycen- 
tcr. The average 3.7-17 keV background count rate in in- 
dividual observations done in 2003-2003 ranges from 7.3 
to 7.8 cts s -1 PCU -1 . Background was not subtracted 
in our timing analyses. 

To determine an appropriate model for spectral anal- 
yses, we used data from both the PCA and HEXTE. 



In the PCA analyses we utilize 2.9-25 keV pulse height 
range data from all layers of PCU 2 only because that 
PCU was operational for all of the observations and be- 
cause its calibration is superior, as judged by the analysis 
of contemporaneous PCA observations of the Crab Neb- 
ula. Spectral extractions and background subtractions 
for both the PCA and HEXTE were performed using 
the "FTOOLS" package, and spectral models were ap- 
plied using "XSPEC" . Both of these software packages 
were provided by NASA/HEASARC. 

3. PULSE TIMING AND ORBITAL ANALYSIS 

A first estimate of the mean pulse period was obtained 
by binning the data from all of the epoch 2 (2002-2003) 
observations in 16-s bins and carrying out an FFT anal- 
ysis of the resultant data train (2 18 points). The FFT 
revealed a highly significant signal at a period of 604.76 
s, as well as 5 very prominent higher harmonics of the 
fundamental. Also evident in the Fourier transform were 
a number of sidebands of the pulsations that were spaced 
in frequency by multiples of 1/4.4 days -1 . These side- 
bands were identified as due to the orbital motion of the 
neutron star about its companion. From the fundamen- 
tal and the prominent harmonics and sidebands, we de- 
rived an average pulse period during these observations 
of 604.689 s. 

We proceeded to fold the data from each of the 71 
pointed observations from 2002-2003 modulo this value 
of the average pulse period. To form a "pulse template" , 
the individual pulse profiles were manually aligned to 
correct for small phase drifts and then averaged (see Fig. 
0). 

We then computed the cross correlation function 
(CCF) between the template and the folded profile from 
each of the individual observations. The phase differ- 
ence for which the CCF reaches a maximum value is the 
best estimate of the pulsar phase relative to the folding 
ephemeris. Pulse arrival time delays are then simply ob- 
tained as the products of the phase differences and the 
folding period. Pulse arrival times were then computed 
by adding each time delay to the time of phase zero of a 
pulse near the middle of the observation interval. This 
analysis yielded a total of 69 pulse arrival times after 
we eliminated the results from 2 of the 71 observations 
for which the CCFs were possibly ambiguous. The pulse 
arrival time delays are listed in Table ^ an d plotted in 
Figure |3 The Doppler delays due to the orbital motion 
are quite apparent in Figure |21 as is an overall quadratic 
behavior due to the slowdown in the rotation rate of the 
neutron star. 

We fit the pulse arrival times with a 7-parameter model 
of the orbit and pulse period behavior. The arrival time 
of the nth pulse is given by 

t n = to + nP H — n 2 PP H — a x sin i cos[f2(i„ — tqq)] 
2 c 

-— a x shiism[2Q,(t n - r 90 ) - w p ] (1) 

where P is the pulse period at time to, P is the pulse pe- 
riod derivative, a^sini is the projected semimajor axis 
of the orbit of the neutron star, r 90 is a reference time, 
which for a circular orbit corresponds to superior con- 
junction, f2 is the orbital frequency which we fix at a 
value of 27r/4.400 radians day -1 , and uj p is the longitude 
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Fig. 1. — PCA counting rate data for the 3.7—17 keV energy band in 8-s time bins for two representative observations of X1908+075, 
i.e., the observations of 2002 December 23 (bottom; Obs. ID 70083-01-01-00) and of 2003 February 8 (top; Obs. ID 70083-02-19-00). The 
vertical lines are drawn every 605 seconds and are roughly aligned with the principal minimum in the pulse profile. The gap in the middle 
of each observation is due to Earth occultation of the source. The non-source background rate is ~ 7.5 cts s _1 PCU . 



of periastron. The last term represents the first order 
term in a Taylor series expansion in the eccentricity e 
and is a reasonable approximation for a mildly eccentric 
orbit. The curve representing the best-fit orbit and pulse 
period behavior is shown superposed on the arrival time 
delays in Fig. [3] The rms scatter of the arrival time de- 
lays about this best fit curve was found to be ~ 6 s and 
was used as an empirical estimate of the uncertainty in 
the individual arrival times. This was the basis on which 
we derived formal confidence limits on the individual fit- 
ted parameters. 

The results of this fit are given in Table [3 The value 
of a x sini of 47.8 lt-s, combined with an orbital period of 
4.4 days, yields a mass function of 6.07 ± 0.35 Mq. This 
points to a fairly massive companion star. The best- 
fit orbital eccentricity is e = 0.021 ± 0.039, and, there- 
fore, the orbit is consistent with being circular. We set 
a 2<T upper limit e < 0.1. The pulse period, after correc- 
tions for orbital motion of the Earth and of the pulsar, 
is determined to be P = 604.684 ± 0.001 s at t =MJD 
52643.3. The pulse period derivative is found to be posi- 
tive, P = (1.22 ±0.09) x 10~ 8 , and thus the neutron star 
rotation is slowing down. Expressed as a fractional rate, 



wc obtain P/P ~ 6.4 x 10~ 4 yr _1 . 

In Fig. 0]we show the measured pulse arrival time de- 
lays for the extensive 2002-2003 epoch 2 data set relative 
to our best-fit model excluding the orbital Doppler shift 
terms. The results are plotted modulo the 4.4-day or- 
bital period. The solid curve is the best fit circular orbit 
with a x smi = 47.8 lt-s. 

A similar pulse phase analysis was then carried out on 
the smaller epoch 1 data set taken approximately two 
years earlier. Only 7 of the pulse arrival time delays 
from this earlier set proved to be useful in the model fit- 
ting, in part because it was difficult to unambiguously 
connect the pulse phase across the large gap between the 
two subsets of epoch 1 observations. The Doppler delays 
for these 7 arrival times, which span only one orbital 
cycle, are plotted in Fig, |S] When fitting these delays, 
we fixed the amplitude (a x sin i) at the value determined 
from the later, more extensive observations, which yield 
much more accurate orbital parameters. The orbital pe- 
riod (P rb) was fixed at the value determined previously 
from 3 years of ASM observations (WRB). The main ob- 
jective in using these observations was to determine a 
relatively accurate orbital phase some two years before 
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TABLE 1 

X1908+075 Pulse Arrival Time Delays 
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70083-02-11-00 


52674.1903 


2560 


9.01 


0.727 


-13.5 


70083-02-12-00 


52674.7445 


1680 


6.08 


0.853 


25.7 


70083-02-12-01 


52674.8198 


1872 


9.37 


0.870 


18.5 


70083-02-13-00 


52675.1768 


2688 


5.23 


0.952 


41.5 


70083-02-14-00 


52675.7355 


2112 


17.71 


0.079 


38.7 


70083-02-14-01 


52675.8025 


2528 


11.68 


0.094 


37.2 


70083-02-15-00 


52676.1636 


2688 


10.81 


0.176 


19.9 


70083-02-16-00 


52676.6994 


1120 


11.96 


0.298 


-18.3 


70083-02-17-00 


52677.0792 


3456 


10.55 


0.384 


-35.3 


70083-02-18-00 


52677.6855 


1232 


10.86 


0.522 


-43.8 


70083-02-19-00 


52678.0684 


3296 


9.68 


0.609 


-33.1 



a Midpoint of observation, MJD = JD - 2, 400, 000.5. 

b Sourcc count rate (2-30 keV) using PCU No. 2. For reference, 1 Crab = 2500 cts s _1 . 

c Orbital phase corresponding to the observation midpoint calculated assuming P or b = 
4.4007 d and the time of phase zero t 90 = MJD 52631.383. 

d Pulse arrival time delay at the Solar System barycenter with respect to a clock with 
constant period P = 604.689 s. 
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Fig. 2. — Pulse profile template used for computing the cross correlation function. The profile is shown for two cycles. The error bar 
shows the approximate uncertainty (±1 a) from variations of the source intensity at a given pulse phase within a typical observation. 
Fluctuations due to counting statistics are negligible. The non-source background rate is ~ 7.5 cts s" 1 PCU -1 . 



TABLE 2 

Results of the Timing Analysis 



Parameter 


2nd-Epoch Parameter a 


lst-Epoch Parameter 1 " 


Comments 


a x sin i (lt-s) 


47.83 ± 0.94 


47.83 (fixed) 




r 90 (MJD) 


52631.383 ±0.013 


51870.06 ±0.15 


superior conjunction for circular orbit 


P(b) 


604.684 ± 0.001 


604.660 ± 0.040 


at MJD 52643.3, 51870, respectively 


P (s s- 1 ) 


(1.22 ±0.09) x 10~ 8 


(fixed) 




e 


0.021 ± 0.039 


(fixed) 


< 0.1 (2a) 


Povb (days) 


4.4007 ± 0.0009 


4.4007 ± 0.0009 


determined from both epochs jointly 


/(M)(M ) 


6.07 ±0.35 


NA 





a Based on 69 pulse arrival times obtained from 2002 December 23 through 2003 February 8. Errors cited 
are single-parameter la confidence limits. 

b Based on 7 pulse arrival times obtained from 2000 November 19 through 2000 November 23. Errors 
cited are single-parameter la confidence limits. 
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Fig. 3. — Pulse arrival time delays for 69 pointed observations made during the interval 2002 December 23 through 2003 February 8. 
The delays are computed relative to a reference pulse ephemeris that has a constant 604.689 s period at the Solar System barycenter. The 
curve shows the arrival time delays predicted by the best-fit model pulsar with a constant pulse period derivative moving in a circular orbit 
(see Tablc|5J. A ±f a error bar, estimated from the rms deviation of the measurements relative to the curve, is plotted in the lower right 
corner. 



the later, more extensive observations. The two orbital 
phase determinations, separated by ~ 2 years, are then 
combined to compute a more accurate determination of 
the orbital period, P orb = 4.4007 ± 0.0009 days. This 
is consistent with, and of comparable precision to, the 
orbital period of 4.400 ± 0.001 days determined from the 
ASM X-ray light curve by WRB. 

By now (2004 April), the ASM has accumulated ~ 
40, 000 individual intensity measurements of X1908+075 
over ~ 8 years. This is more than twice the amount of 
data contained in the original light curve of WRB. We 
have therefore redone the light curve analysis using all 
the presently available ASM data to derive a new X-ray 
intensity-based value for the orbital period. Our result is 
P mh = 4.4006 ± 0.0006 days. This is consistent with the 
earlier result of WRB and with our pulse-timing-based 
result as well. The folded ASM light curves are shown in 
Fig. El 

4. ORBITAL PHASE-DEPENDENT SPECTRAL ANALYSIS 

The X-ray intensity of X1908+075 is strongly modu- 
lated with orbital phase. A qualitative feeling for the 
energy-dependence of the modulation can be obtained 



from the folded ASM data (Fig. EJ). The ASM light 
curves are extremely well sampled both as a function of 
orbital phase and over many orbital cycles. We define 
the degree of modulation as (I max - I mi „)/(/ mM + I m in), 
where 7 max and I m m are the maximum and minimum 
intensities observed in the light curves. The degrees of 
modulation are found to be 0.41 ± 0.18, 0.70 ± 0.15, and 
0.38 ± 0.04 for the 1.5 - 3, 3 - 5, and 5 - 12 keV en- 
ergy bands, respectively. The uncertainties have been 
estimated primarily on the basis of statistical fluctua- 
tions in 7 m i n , and do not include any contributions due 
to possible systematic baseline offsets. Nonetheless, we 
take these numbers, as well as the different shapes of the 
folded light curves, to be evidence of energy dependence 
of the orbital modulation. 

In spite of the fact that the orbital light curves pro- 
duced with the ASM are very well sampled, the energy 
resolution (limited to 3 channels) and calibration are not 
good for determining detailed spectral parameters. The 
PCA data are better suited to this purpose. 

To help guide the spectral analysis, we consider the 
overall average spectrum of X1908+075 during the epoch 
1 and epoch 2 RXTE pointed observations. In each case, 
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Fig. 4. — Doppler delay curve for the same set of observations as shown in Fig. [3] but plotted as a function of orbital phase. Each 
measurement is plotted twice in order to show two complete orbital cycles. The quadratic term representing the change in intrinsic pulse 
period, evident in Fig. |!|] has been removed. The curve represents the best-fit circular orbit. 



the PCA spectrum (2-25 keV; PCU 2) shows evidence of 
low-energy attenuation by a substantial column density 
and the presence of line emission due to neutral or ionized 
Fe at energies of approximately 6.4-6.7 keV. The HEXTE 
spectrum (both clusters) in the range 15-70 keV shows 
a high-energy cutoff that rules out a simple power-law 
model. For both epochs, thermal bremsstrahlung and 
cutoff power-law models yield reasonable fits of the con- 
tinuum spectrum. We proceed with the bremsstrahlung 
model since it has fewer free parameters, gives results 
for each epoch that are remarkably consistent, and (in 
hindsight) gives best-fit temperatures from the individ- 
ual observations that agree well with each other. Since 
the spectral resolution of the PCA is at best marginal 
for determination of the Fe line energy, we assume that 
any line emission near these energies is at 6.7 keV; this 
should be an adequate model of the Fe line region for 
the present purposes. The fit of the average spectrum 
from epoch 1 (35 ks of exposure) yields a column den- 
sity N H = (1.12 ± 0.02) x 10 23 cm" 2 , a temperature 
kT = 23.3 ±0.5 keV, and a value of chi-square per degree 
of freedom \t — 2.42. The fit of the spectrum from 153 
ks of exposure in epoch 2 yields N H = (1.22±0.02) x 10 23 
cm" 2 , kT = 22.6 ±0.3 keV, and xl = 4.84. These values 



for xt are formally unacceptable. However, the average 
spectra comprised data obtained from large accumulated 
exposure times so that the uncertainties from counting 
statistics are very small. In fact, in our analysis of the 
PCA spectral data we estimated the error for each en- 
ergy bin by adding the statistical error in quadrature 
with a presumed 1% systematic error which represents 
the uncertainty in th e instrument response function (e.g., 
ISobczak et alJl2000|) . For much of the energy range these 
presumed systematic errors dominate the statistical er- 
rors. Perhaps most importantly with regard to the x 2 
values of the fits of the average spectra, the Njj values 
vary significantly from observation to observation. In 
such model based on a single value for Njj can- 

not be expected to precisely fit the average spectrum. 

Next, we analyzed the net spectrum for each of the 69 
epoch 2 pointed observations for which we have deter- 
mined a pulse arrival time. The analysis was repeated 
for each of three spectral models in order to estimate how 
the formal statistical uncertainties compare with those 
due to the particular choice of model. Each model is in- 
tended to provide a simple spectral shape that is likely to 
fit the data well and to yield a good estimate of the low- 
energy absorption. All the models involved a thermal 



LEVINE et al. 



60 




-60 r , , , i , , , i , , , i , , , i , , , i 

0.0 0.2 0.4 0.6 0.8 1.0 

Orbital Phase 



Fig. 5. — Seven pulse arrival time delays for observations during 2000 November (approximately two years prior to the more extensive 
observations shown in Fig. EJ. These data spanned only one orbital cycle of X1908+075. The curve represents the best-fit circular orbit 
(see text). The uncertainty in the arrival time delay (±1 <r), estimated from the rms deviation of the measurements relative to the curve, 
is shown as an error bar on one of the measurements. 



bremsstrahlung spectrum, a 6.7 keV iron line (which, as 
noted above, also represents line emission near 6.4 keV), 
and absorption by the interstellar medium (ISM) and a 
stellar wind. 

Inspection of the raw pulse height spectra reveals that 
there are 'excess' counts at low energy (i.e., < 6 keV) 
with respect to the energy dependence expected from 
photoelectric absorption in a neutral gas. Two possible 
explanations of this form of the low-energy portion of the 
spectra immediately arise. First, the observed spectrum 
could be formed by both X-rays that propagate directly 
to the detectors along paths that pass close to the com- 
panion star and therefore are heavily attenuated, as well 
as X-rays that are scattered farther out in the stellar 
wind and are not attenuated so heavily by photoelectric 
absorption. Second, if the stellar wind is partially ion- 
ized, then the energy dependence of any absorption will 
not be as strong as that for strictly neutral gas. The 
combination of these effects can be approximately taken 
into account by allowing the presence of two spectral 
components that are attenuated by different amounts. 

We take the intrinsic emission spectrum of the source 



to be 

S(E) = dBr(E, kT) + C 2 G{E) (2) 

where Br(E, kT) represents an optically thin thermal 
bremsstrahlung spectrum characterized by temperature 
T and G(E) represents a Gaussian profile line centered 
at 6.7 keV with width AE(la) = 0.3 keV. This is used, 
in turn, in a generic spectral model: 

F x =ex V [-N HiaM a(E)] x 

{C 3 exp[- N Hwiad a(E)] + (1 - C 3 )}S(E) . (3) 

C\ , Ci , and C3 are parameters of the model. Each of the 
three spectral models that we apply to the observations 
uses a particular form of this generic model. In Model 1, 
the thermal bremsstrahlung spectrum is attenuated by 
the ISM and a stellar wind, kT is free to vary from one 
observation to another, and C3 is assigned a value of one. 
The ISM and wind column densities are not determined 
separately; only their sum is determined. Model 2 is the 
same as Model 1, except that kT is fixed at the value of 
22.6 keV that was determined by the fit of the average 
spectrum of the epoch 2 observations. Model 3 is similar 
to Model 2, except that the absorption is characterized 
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Fig. 6. — Average intensity as a function of orbital phase pro- 
duced by folding ~ 8 years of ASM results at the period P = 4.4006 
days. The dotted, solid, and dashed histograms correspond to the 
(nominal) 1-3.5, 3.5-5, and 5-12 keV energy bands. Each curve is 
shown for two cycles. For each of these energy bands, an intensity 
of 1 ASM count s — 1 corresponds to ~ 40 mCrab. 



by 'partial covering' and so the parameter C3 is deter- 
mined for each observation. For Models 1 and 2, the free 
parameters were simply determined by varying the free 
parameters so as to minimize the x 2 statistic. For Model 
3 such a procedure did not yield acceptable results, so an 
iterative procedure was used. In this procedure, the value 
of Njj lsu was fixed and the remaining free parameters in- 
cluding -/Vff w were determined for each observation by 
a minimum x fit with XSPEC. A fit of a model of the 
orbital phase dependence of iVij wind was then made to 
the resulting NH wiad values to determine wind parame- 
ters (with (3 = 0, see §5 below). Then the spectral fit of 
each observation was redone using XSPEC with Nh ism 
as a free parameter and with the value of Nn vind fixed to 
the value predicted from the fit of the wind parameters. 
The weighted mean of the Nh 1sm values then was used as 
a revised estimate of this quantity. This entire procedure 
was iterated 4 times to obtain convergence. In the final 
iteration, Nh 1sm was set to 4.6 x 10 22 atoms cm~ 2 , which 
indicates that substantial absorption occurs in the ISM 
(although it should be noted that this column density is, 
at all orbital phases, comparable or smaller than that in 
the stellar wind). 

The statistics of the individual observations were not 
sufficient to conclude that one of the spectral models 
fit the data significantly better than the other models. 
Rather, all three models yield comparable quality fits 
with reduced x 2 values in the range 0.5-2.0 for all but 
a few observations. As might be expected, the value of 
the normalization constant C\ does not appear to be de- 
pendent on orbital phase, but the iron line intensity does 
appear to depend on orbital phase by as much as 30% in 
the Model 2 results and less in the Model 1 and Model 
3 results. In the Model 3 fit results, the value of C3, the 



covering fraction, is always in the range 0.4-1.0 and, for 
~ 85% of the observations, is in the range 0.85-1.0. 

The total neutral hydrogen column densities (Nh = 
-^ffisM + -^ffwmd ) f° r the 69 observations are shown as a 
function of time in Fig. the top panel shows results 
from spectral Model 1, while the bottom panel shows 
those for the partial covering model (Model 3). Note 
the pronounced variation of Nh with maxima at the ex- 
pected orbital phases. A comparison of the two panels 
provides an idea of the sensitivity of Nh to the choice 
of spectral model; spectral Model 2 yields qualitatively 
similar results. The fitted values of Nh based on Mod- 
els 1 and 3 are also plotted modulo the orbital phase in 
Figure Q3 

The average spectrum from seven observations per- 
formed at orbital phases when the column density is 
expected to be high, i.e., \(j> or b\ < 0.12, is shown in 
the top and middle panels of Figure |5| The degree of 
low energy attenuation contrasts with that of the spec- 
trum from eight observations performed at orbital phases 
when the column density is expected to be low (Fig. [5] 
lower panel). We have fit each average spectrum with 
both Models 2 and 3, but we fixed the value of the 
bremsstrahlung temperature in each case to that deter- 
mined from the average spectrum of all the epoch 2 ob- 
servations. We also fixed the energy of the Gaussian line 
center at 6.7 keV and the line width at au ne = 0.3 keV. 
The results are given in Table [3] The Model 2 best- 
fit high Nh spectrum and the Model 3 best-fit low and 
high Nh spectra are also shown in Fig. ED The plots 
illustrate the low energy attenuation and its variation. 
Spectral features are quite apparent at ~ 7 keV; these 
comprise Fe K line emission and absorption edges. The 
iron absorption is particularly prominent in the high Nh 
spectrum. Both the low and high Nh spectra are fit 
somewhat better by Model 3, although neither model fits 
the data well; the values of reduced x 2 are significantly 
above 1. 

From the spectral fits, we find that the X-ray flux, av- 
eraged over an observation and adjusted for low-energy 
absorption, is most often in the range Fx (2-30 keV) 
~ 3-8 xl0~ 10 ergs cm~ 2 s _1 . Thus the X-ray lumi- 
nosity of X1908+075 is typically L x (2-30 keV) ~ 2-6 
xlO 36 (D/8 kpc) 2 ergs s" 1 where D is the distance to 
the source, and, even at the high end of this range, is 
far below the Eddington limit for a 1.4 Mq neutron star, 
i.e., Lndd = 2.0 x 10 38 ergs s _1 . This is consistent with 
the picture that the neutron star is accreting from a wind 
from the companion star rather than being fed via Roche 
lobe overflow. 

5. ABSORPTION IN A STELLAR WIND AND THE BINARY 

SYSTEM 

The variation of the column density with orbital phase 
is caused by the movement of the X-ray source through 
the stellar wind of the companion star which is likely 
a spatially and temporally complex medium. In order 
to estimate the orbital inclination as well as other or- 
bital parameters and properties of the companion star, 
we proceed by modelling the column density variation 
as orbital-phase dependent absorption in a spherically 
symmetric and temporally constant wind from the com- 
panion star. 

The wind from an early-type star is often mod- 
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Fig. 7. — Estimates of Njj = Ng ind + Nn 1SM derived from spectral fits plotted as a function of time expressed in units of orbital 
cycles. Integral values of time, i.e., tic' times of orbital phase zero, correspond to superior conjunction of the neutron star. The upper 
plot shows estimates derived from fits of spectral Model 1 with uncertainties (±1 <x) taken from the XSPEC fits. The lower plot shows 
estimates from spectral Model 3. with uncertainties (±1 <x) taken from those for Nh ind from the final set of XSPEC runs since we assume 
the uncertainty in Njj ism may be neglected in this case. The curves show the best fits of the orbital-phase dependent column density with 
j3 = 0; the best-fit curve for j3 = 1 (see text) is nearly indistinguishable from the curve shown here. For this plot, the long gap between 
sets of observations was artificially shortened by 4.0 orbital cycles. 



elled as a steady-state, spherically symmetric flow at 
a velocity that gradually increases from zero to a 
terminal velocity of order of magnitude ~ 1000 km 
s _1 . The wind density profile then follows from 
the assumption of a radius-independent mass flux 
lILucv fc SolmnT^ll97fl |Castor. Abbott, fc Kleinlll97fit 
ILamers fe Cassinellilll99SA and references therein'! . For 
our purposes the radial flow velocity may be taken to be 



v(r) 



■Xl-Rc/r) 



13 



(4) 



where r is the distance from the center of the companion 
star, Voa is the wind terminal velocity, R c is the radius of 
the companion star, and is often in the range 0.7-1.2 
for early-type st ars fe.g.. iGroenewegen fc Lamerslll989t 
IPuls et al.lll996fl . The wind density profile then follows 
from the conservation of mass: 

n (l - R c /a) 



13 



n(r) 



(5) 



(r/a) 2 (l - R c /r)P 

where ii{r) is the number density of hydrogen atoms in 
the stellar wind at the distance r from the companion 



star and no is the number density at the distance a. For a 
circular orbit with orbital inclination angle i, the instan- 
taneous column density of material between the neutron 
star and the observer, Nh , is given by: 



N 



H : 



N, 



n[r(s)]ds (6) 



■■N Hls 



an (l - R c /a) 



ds' 



4J) 



io r' 2 (l- R'Jr'Y y 

where Nh 1sm and Nn^ iad are the separate contributions 
from the interstellar medium and the stellar wind, re- 
spectively, s is the distance along the line from the neu- 
tron star toward the observer, and primed quantities are 
normalized relative to a. 

The radial distance r is a function of the distance s 
and the angle (f> subtended at the neutron star between 
the direction to the observer and the radial direction. In 
terms of normalized quantities, we have 



1 + s' 2 - 2s' cos <f> . (8) 
in turn, is related to the inclination angle 
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Fig. 8. — (top) Njj estimates for spectral Model 1 as a function of orbital phase. The points shown are those shown in the top panel of 
Fig. 17] ( bottom) Njj estimates for spectral Model 3. In both panels, the solid curve shows the best-fit /3 = stellar wind model. As noted 
in Figwl the best-fit curve for the /3 = 1 model is nearly indistinguishable. The points and curves are shown in each of two cycles. 



TABLE 3 
Spectral Fit Results 





Low Njj 


Low Njj 


High Njj 


High Njj 


Parameter 


Model 2 


Model 3 


Model 2 


Model 3 


Ci a 


0.1012 ± 0.0004 


0.1023 ±0.0005 


0.0880 ± 0.0004 


0.0905 ± 0.0005 


C 2 b 


5.7 ±1.1 


6.1 ± 1.1 


3.3 ± 1.0 


4.8 ± 1.1 


^ism C 




4.7 (fixed) 




4.7 (fixed) 


Njj ■ a c 




8.5 ± 1.3 




24.7 ±0.7 


AT H ^total) c 


9.3 ±0.1 




24.3 ±0.2 




Cz (covering fraction) 




0.66 ±0.06 




0.916 ± 0.008 


x 2 


126.2 


115.0 


222.3 


142.5 



Degrees of freedom 

xl 



46 
2.74 



45 
2.56 



16 
4.83 



45 
3.17 



Note. — All quantities are defined in eqs. 2 and 3. For all 4 fits, the bremsstrahlung 
temperature was fixed at kT = 22.6 keV, the spectral line center energy was fixed at 6.7 keV and 
the line width fixed at a = 0.3 keV. The low Njj spectrum is the average from 8 observations. 
The high Njj spectrum is the average from 7 observations (see text). 

formalization factor for XSPEC component "bremss" to give flux in units of photons cm -2 
s- 1 keV" 1 

b Iron line flux in units of 10 — 4 photons cm -2 s _1 
c In units of 10 22 H-atoms cm -2 



and the orbital phase by: 

cos 4> = — sini cos[f2(i — rgo)] 



(9) 



For the simple case of an inverse square law density pro- 
file, i.e., /3 = 0, the wind column density integral equals 



12 



LEVINE et al. 





-2.0 






1 1 1 


men in |_| 




-2.5 








Model 2 : 




-3.0 


- -P 










-3.5 










> 

o 














-2.0 








High N H 


CO 


-2.5 








Model 3 : 


E 

o 


-3.0 










Photons 


-3.5 










cn 


-2.0 








Low N H : 


o 

_i 


-2.5 

-3.0 

-3.5 
-4.0 








Model 3 : 




3 




5 


10 


20 



Energy (keV) 



Fig. 9. — Inferred photon number spectra (points with error bars) 
of the accumulated data from 7 observations when the absorbing 
column density was inferred to be relatively high (top and middle 
panels), and of the data from 8 observations when the column 
density was relatively low (bottom panel). The histogram in the 
top panel shows the best-fit of the Model 2 spectrum while the 
histograms in the middle and lower panels show the best fit of the 
Model 3 spectrum. 



ano4>/ sin^> if (f> is given in radians. In this case, the wind 
column densities, in units of noa, are 1, tt/2, 371-72/4, and 
oo for = 0, 7r/2, 37r/4, and it radians respectively, and 
the ratio of maximum to minimum value is given simply 
by ^(maxVATj^H = (tt/2 + i)/{jz/2 - %). 

Since we know tqq from our orbital analysis, there are 5 
parameters that can potentially be determined from fits 
to the orbital modulation of Nh- i, (3, R c , Nh 1sm , and 
the product an^(\ — R c /a)". In the analysis procedure, 
for each of a set of assumed trial values for i, (3, and R c , 
we performed a linear least squares fit of the model de- 
fined by eq. (7) to the column densities determined from 
our Model 1 spectral fits, and thereby obtained values for 
Nh ism and ano(l — Rc/a) 13 as well as an estimate of the 
root-mean-square deviation of the column densities from 
the fit. We have fit models for a wide range of orbital 
inclinations. For each trial value of i, we use our best-fit 
value for f(M) (see Table|2J) and assume that the mass of 
the neutron star is the canonical IAMq to compute the 
mass of the companion M c . Using the resulting mass ra- 
tio, i, and our best fit value for a x sin i, we then compute 
the orbital separation a. For the given inclination, the 
set of trial values of R c covers a range up to a maximum 
value taken to be the smaller of the Roche lobe radius of 



the companion star or the radius which would produce 
a grazing eclipse. Finally, for simplicity, we considered 
only three discrete values for the wind parameter (3: 0, 
1/2, and 1. 

Since we do not have any reliable way to estimate the 
uncertainties in the measured column densities due to 
systematic errors, we have taken the rms spread between 
the values found from the spectral analysis (for the 69 
individual observations) and the best-fit model to be the 
Ic7 error. When we do this, the minimum value of x 2 
should be 65, by definition, since there are 69 data points 
and 4 fitted parameters (not counting 0), The contour at 
X 2 ~ 65 + 4.6 = 69.6 is then relevant for estimation of the 
90% confidence region in the i — R c plane (with Nh ism 
and ano{l — R c /a)^ treated as uninteresting parameters). 
In Fig. 1101 we show contours of x 2 m this plane from 
fits of the (3 = 1 model. We note that, in the limit 
of vanishing R Cl the value of (3 becomes irrelevant and 
the model is equivalent to the j3 = model. Thus, the 
contours at small R c in Fig. 1101 indicate the acceptable 
range of inclinations in the /3 = fits. We find that the 
best fits for each of the three values of (3 are very similar 
in quality, so we cannot discriminate among these values 
of (3 using this criterion alone. 

The best fits of the (3 = 1 model to the values of Nh 
from spectral Models 1, 2, and 3 give rms observed- 
minus-calculated values of 4.1 x 10 22 , 5.1 x 10 22 , and 
7.7 x 10 22 H atoms cm -2 , respectively. The (3 = 
and 1/2 cases give slightly larger residuals. The best 
fit to the spectral Model 1 results of a (3 = model is 
characterized by i = 64.5°, M c = 11M Q , a = 60 lt-s, 
N H 1SM = 4.1 x 10 22 cm" 2 , and N Hvind (<f> = 0) = an Q 
-1.7 x 10 22 cm -2 . The best fit to the spectral Model 1 
results of a (3 = 1 model has i = 48.5°, M c = 17M©, 
R c = 16i? Q , a = 69 lt-s, N HlSM = 5.0 x 10 22 cm" 2 , and 
iVff wtad (0 = 0) = 3.6 x 10 22 cm~ 2 . The column densities 
predicted by these models are practically indistinguish- 
able from each other; they are shown in Figs. [7| and 
|S| We emphasize that acceptable fits are obtained for a 
considerable range of parameters, and that the best-fit 
parameters given above are merely illustrative of the al- 
lowed values. Our estimates of the system parameters 
for X1908+075 are summarized in Tabled 

The results of the column density fits are shown as con- 
straints on the companion star mass and radius in Fig. 
im and allow us to set, e.g., assuming (3 = 1.0, upper 
limits of ~ 30M Q and ~ 22i? Q on the mass and radius, 
respectively, of the companion star. We note that for 
M c > 13.5-Mq, the upper limits on the companion ra- 
dius correspond to the Roche lobe radius, while the lack 
of an apparent X-ray eclipse bounds the allowed radii for 
M c < 13.5M Q . In Fig. ^2 we also show the measured 
masses and radii of the companion stars in 6 well known 
high-mass X-ray binaries and mass-radius relations for 
zero-age main-sequence (ZAMS) stars, for stars that have 
exhausted half of the hydrogen at their centers, and for 
stars at the end of the main sequence phase of their evolu- 
tion, i.e., the terminal-age main sequence (TAMS), where 
the mass is that indicated on the axis (Ph. Podsiadlowski 
2004, private communication). These mass-radius rela- 
tions were derived for single (isolated) stars neglecting 
the effects of any mass loss. 

Each set of system parameters and wind density pro- 
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Fig. 10. — Contours (solid curves) of \ m the i — R c plane from fits of the column densities for wind density profiles with /3 = 1 to the 
spectral Model 1 Njj measurements shown in the upper panel of Fig. [7] The three contours show the values for which \ 2 = (1 + A / ^65)x^j n 
where A = 0.25, 4.0, or 4.61. The outer contour corresponds to the formal 90% confidence limit for two interesting parameters. For each 
value of the inclination, the mass of the companion star is determined to within a small uncertainty. The long-dash curve represents the 
radius of the Roche lobe of the companion star, the medium-dash curve represents the radii where grazing eclipses would occur, and the 
three short-dash curves represent the radii of stars on the zero-age main sequence (ZAMS), stars that have exhausted half of the hydrogen at 
their centers, and stars on the the terminal age main sequence (TAMS; see text). The mass-radius relations were provided by Podsiadlowski 
(2004, private communication). 



file that is consistent with the column density measure- 
ments yields a wind mass flux if we assume a wind 
terminal velocity. We u s e the prescription given in 
IVink. de Koter. fc Lamersl lj2000|) . which predicts wind 
terminal velocities on the basis of the stellar mass, stel- 
lar radius, effective temperature, a nd luminosity (see also 
lLamers. Snow, fc LindnolmH l995). In turn, we interpo- 
late the stellar evolution calculations of Podsiadlowski for 
main sequence stars to estimate effective temperatures. 
Thus, we have not carried out this calculation for stellar 
radii smaller than those on the zero-age main sequence or 
larger than those on the terminal-age main sequence (see 
Fig. 1111) . We can compare these "observed" wind mass 
fluxes with those expected for a star with given mass, 
effective temperature, and luminosity. The "expected" 
mass flux is computed using formulae also presented by 
Vink et al. These formulae are based on empirical fits 
of mass fluxes inferred from ultraviolet, radio, and Ha 
observations of many O and B stars. The same wind ter- 
minal velocity which we use to obtain an estimate of the 



observed mass flux is used in the Vink et al. prescription 
for wind mass flux. The mass fluxes used by Vink et 
al. in the empirical fits are for single stars, and are not, 
in particular, for stars in binaries which nearly fill their 
Roche lobes. 

We find that the "observed" wind mass loss rate is 
larger in all cases than the "expected" mass loss rate. In 
Fig. 1111 medium-size dots show those combinations of 
system parameters in which the "observed" wind mass 
loss rate exceeds the "expected" rate by a factor of 3 or 
less. In the other cases for which we computed the "ob- 
served" wind mass loss rate, it exceeds the "expected" 
rate by more than a factor of 3, and often by more than 
a factor of 10. Most of these dots apply only to the (3 = 1 
case; one dot applies to the (3 = 1/2 case; none apply to 
the (3 = case. 

The "expected" mass flux most closely approaches the 
"observed" mass flux for the case with (3=1, M c = 
22.8M Q , R c = 20.0i? Q , ~ 800 km s^ 1 for which 
we obtain M bs ~ 2.1 x 10~ 6 Af Q yr -1 and M exp ~ 1.5 x 
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TABLE 4 

System Parameters 



Parameter 



P = = 1/2 



= 1 



a (lt-s) 
M c (M Q ) 
R c (Rq) 

N H ^. nd (<j> = 0) a (10 22 cm" 2 ) 
M (10 -6 Af Q yr" 1 ) 



54° - 72° 43° - 70° 38° - 68° 

58 - 65 58 - 75 59 - 81 

9 - 14 9 - 22 10 - 31 
4 - 15 4 - 19 4 - 22 

3.4 - 6.8 3.1 - 6.9 2.5 - 6.8 

2.3 - 15 1.9 - 15 1.3 - 14 



Note. — Each parameter range is denned by x 2 < (1 + 
4.61/65)x 2 ra j„ for the fits of eq. 7 to the measured column densi- 
ties. The parameter ranges are also restricted by the requirements 

that R c < Reclipse, Rc < RrocHc, an d by RzAMS < Rc < RtAMS 

for the associated value of M c , i.e., the radius of the companion star 
is required to be in the range of rad ii o f th e m ain sequence stars 
calculated by Podsiadlowski (see Figs. I1UI and llll . 

a The column density of the wind in the radial direction outward to 
infinity from the position of the neutron star. 



10 _6 Mq yr _1 . Note that this case is not realistic because 
the companion radius is equal to the radius of its Roche 
lobe. It is likely that the companion underfills its Roche 
lobe; otherwise the mass accretion rate would be higher 
and the X-ray luminosity would be expected to be at the 
Eddington limit. 

6. DISCUSSION 

We have found 605 s pulsations in the X-ray intensity 
of X1908+075. Doppler shifts of the pulse frequency, and 
changes in the intensity and low-energy attenuation that 
cyclically recur at the previously known orbital period of 
4.4 days allow us to measure orbital parameters and to 
conclude that the system contains a highly magnetized 
neutron star orbiting in the wind of a massive companion 
star. 

The observed variations in the shape of the spectrum 
and in the overall intensity indicate that the X radia- 
tion from the neutron star is both absorbed and scat- 
tered in the wind. We have estimated the absorption 
on the basis of a crude spectral model that assumes 
that any attenuation is caused by photoelectric absorp- 
tion in neutral (unionized) gas with solar element abun- 
dances. In this model, we did not include the effects 
of the photoionization of the wind by the X-rays. We 
can now use our estimates of the wind column density 
and the X-ray luminosity to roughly estimate the ioniza- 
tion p arameter £ = Lx/nr 2 l|Tarter. Tucker, fe Sakoeterl 
1969) at a typical location in the binary system, where 
n is the number density of atoms in the wind and r is 
the distance from the X-ray source. If we take the or- 
bital radius a as a characteristic distance within the bi- 
nary system, i.e., we set r — a, and for a density use 
n ~ N Hvriiid ((/) = 0)/ a, we obtain £ = L x /(N H a) ~ 20- 
60 ergs cm s . Under these conditions we would expect 
that most atoms would lose their outer electrons, but 
atoms of carbon, oxygen, etc. would retain at least their 
K-shell electrons llTarter et alll969HKallman fc McCravl 
11982ft . iKrolik fc Kallmanl lj!98^T have estimated the 
equilibrium conditions of gases photoionized by contin- 
uum X-ray spectra. They consider the condition of the 
gas as a function of a modified ionization parameter 
S = Lx I (4ircnr 2 kT) where c is the speed of light, k 



is Boltzmann's constant and T is the gas temperature. 
For a characteristic location in the X1908+075 system, 
we have S ~ 8(T/10 5 K)~ 1 . Fo r each of the two spectral 
shapes that IKrolik fc Kallmanl (^984) consider in detail, 
one can roughly sketch the parameters of a self-consistent 
solution taking into account that the actual density will 
be higher than that estimated on the basis of the neu- 
tral matter absorption cross section. The parameters 
of the solutions in the two cases are roughly the same: 
T ~ 4 x 10 4 K, 5 ~ 10, the X-ray opacity in the ~2-6 
keV band is reduced by a factor of ~ 2 from the value 
expected for cold neutral gas of cosmic abundances, and 
A# wind (0 = 0) ~ 1 X 10 23 cm~ 2 . Since, in this regime, 
the temperature is a steep function of 2, it is likely to 
vary over a wide range at different places in the system. 
A more realistic estimate of the physical properties of the 
wind would need to self-consistently consider the actual 
source X-ray spectrum, the wind velocity, the possibility 
of non-solar abundances, the position-dependent degree 
of ionization in the wind, the effects of the ionization 
upon the wind acceleration, and the transfer of radiation 
in the system including both the effects of scattering and 
absorption. This is beyond the scope of this paper. 

Using our fitted column de nsities and wind speeds es- 
timate d from the formula of iVink. de Koter. fc Lamersl 
(2000), we find the mass loss rate in the wind from the 
companion star to be > 1.3 x 10 _6 M Q yr _1 for all the 
cases with acceptable values of x 2 an d with companion 
star radii in the range expected for main sequence stars. 
For those cases where the estimated mass loss rate is not 
more than a factor of order 3 greater than an empiri- 
cal prediction based on the stellar mass, temperature, 
and luminosity, we find the rate must be > 2 x 10~ 6 M Q 
yr _1 . The depression of the X-ray opacity of the gas from 
ionization implies that the wind mass loss rate must be 
higher than this latter rate, likely, as discussed above, by 
a factor of two and possibly by a larger factor. 

As noted in §5, over the allowed region of the M c — R c 
plane that we have identified, such mass loss rates 
are significantly higher than those predicted using the 
empirically-based mass loss rate prescriptions of Vink et 
al. However, one should note that these prescriptions do 
not take into account any of the effects of the star being 
in a binary system, including the presence of a critical po- 
tential lobe. Furthermore, we used the effective temper- 
atures computed by Podsiadlowski for stars that evolve 
without mass loss and without being affected by any of 
the phenomena that occur in a binary system. Thus, we 
cannot exclude the possibility that the companion star 
is on the main sequence. 

Another possibility is that the companion is a Wolf- 
Rayet (WR) star. A WR star could well have a mass 
that is consistent with the value that we find for the 
companion, but its radius would be much smaller than 
that of a main-sequence star o f comparable mass. It 
would have a prodigious wind l)Nugis fc Lamera 2000) 
that could reach or surpass the mass loss rate that we 
infer (with assumption s) for X1908+075. According to 
iNugis fc Lanie rs (2000), WR stars of either type WN 
or WC and mass in the range ~ 9-15 M Q have wind 
mass loss rates between 4 and 20 x 10 -6 M Q yr" 1 . Of 
course, the winds of WR stars, which have little or no hy- 
drogen, would have significantly larger X-ray absorption 
cross sections per unit mass than winds from unevolved 
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Fig. 11. — Limits on the companion star radius as a function of the companion mass determined from fits of the variation of column 
density with orbital phase. The solid, short-dashed, and l ong -dashed lines show limits from the (3 =0, 1/2, and 1 fits, respectively, and 
correspond to the formal 90% confidence contours (cf. Fig. 1101 . High values of M c correspond to low values of i and vice versa. The lower 
bounds on the radius are arbitrarily set at ~ 1 Rq . The medium-size dots show those fits fo r which the wind mass loss rate in ferred from 
the observations is no more than a factor of 3 greater than the rate predicted by the recipe in Vink. dc Kotcr, & Lamers (2000), (see text). 
The dotted curves are mass-radius relations for main-sequence stars where the mass is that indicated on the horizontal axis; the lowest 
curve represents the zero-age main sequence, the middle curve represents stars that have exhausted half of the hydrogen at their centers, 
and the upper curve represents the terminal age main sequence (see text). The mass-radius relations were provided by Podsiadlowski (2004, 
private commu nication). Also shown on the plot are the measured masses and radii of the "normal" stars in six well known high-mass 
X-ray binaries (Joss & Rappaport 1984; Nagasc 1989; A. Levine 1984, private communication). The associated error bars crudely represent 
90% confidence uncertainties. 



stars. If there is little hydrogen in the wind, we may 
have overestimated the wind mass loss rate, but such a 
composition would nonetheless clearly indicate the WR 
nature of the companion. 

If the companion star is indeed a WR star, then the 
X1908+075 system has potentially important implica- 
tions, in general, for binary stellar evolution, and for 
the formation, in particular, of neutron star-black hole 
(NS-BH) systems. If it is a WR star, then it is likely 
the He or CO remnant core of a star that was originally 
much more massive, possib ly of mass ~ 23 — 35 Mq (e.g., 
iHurlev. Pols, fc ToutJl2000j) . In this case we expect the 
companion to undergo core collapse in 10 4 to 10 5 years 
and leave behind a stellar-mass black hole ijBrown et al.l 
2001). Thus, there is a possibility that the X1908+075 
system is the progenitor of a neutron star-black hole 
binary where the neutron star formed first. This would 
be significant in at least two ways. First, in spite of 



the fact that six NS-NS binaries have been discov ered 
(see, e.g.. lBurgav et af]l2003l: TChampion et alJfeOO^ . no 
NS-BH binaries ha ve yet been found. There are theo- 
retica l arguments iSipior. Portegies Zwart fc Nelemansl 
12001 iPfahl. Podsiadlowski. fe RappaportJ l2004|) which 
suggest that the latter binaries should be at least a 
factor of 10 less populous th an their NS-NS cousins 
(but see iVoss fc Taurisl l2003j) . An identified progen- 
itor would help theorists better estimate the current 
population of NS-BH binaries. Second, if NS-BH 
binaries are formed, the issue of which collapsed star 
forms first is also important to our understanding 
of binary stellar evolution. The most direct way 
of producing these systems is for the BH to form 
first (from the mo re massive component of the sys- 
tem; see, e.g.. iPortegies Zwart. Verbunt. fc Ergma 
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2001 ILee. Brown, fc Wiiersl 12001 

Podsiadlowski. Rappaport. fc Hani l2002fl . but there 
are also channel s where the NS form s first, as discussed 
above (see, e.g.. lVoss fc Taurisll2003j) . 

Since this p ap er was originally written, 
iMorel fc Orosdidierl ((2004) have reported the re- 
sults of near-infrared observations of stars in or close 
to the error box of X1908+075, and have found a star 
whose JHK magnitudes and colors, and H and K band 
spectra, suggest an O or B supergiant at d ~ 7 kpc. 
They believe that this star is likely to be the counterpart 
of X1908+075. If confirmed, this identification would 
rule against a WR companion star. We expect to obtain 



observations of X1908+075 with the Chandra X-ray 
Observatory, and to thereby obtain an X-ray position 
sufficiently accurate to secure the optical identification. 
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